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Abstract. In this paper, we study a distribution T of labeled particles 
{1, 2,...,n} on a continuous ring. It arises in three different ways, all 
related to the multi-type TASEP on a ring. We prove formulas for 
the probability density function for some (cylic) permutations and give 
conjectures for a larger class. We give a complete conjecture for the 
probability of two particles i,j being next to each other on the cycle, 
for which we prove some cases. 

We end with observations of the similarities between the TASEP on 
a ring and the well-studied Razumov-Stroganov process. 


1. Introduction 

In this paper, we study a distribution T of labeled particles on a contin¬ 
uous ring. We call this distribution continuous TASEP on a ring, and it 
arises in three different ways: 

(1) As the limit of the stationary distribution of totally assymetric ex¬ 
clusion process (TASEP) on a ring. 

(2) As the projection to the last row of a random continuous multiline 
queue. 

(3) As the stationary distribution of the so-called (continuous) process 
of the last row. 

Exact definitions are given in Section [2j The equivalence of these three 
descriptions already in the discrete case follows from the seminal work con¬ 
ducted by Ferrari and Martin [12] . The first two explicitly and the third 
implicitly, as described in Section 12.31 The limit of the TASEP considered 
here keeps the number of jumping particles constant and let the number of 
vacant positions tends to infinity. 

In previous works of several authors e.g. the finite 

version of this TASEP has been proved to have many remarkable properties 
and to be connected to other mathematical objects: shape of random n -core 
partitions, random walks in an affine Weyl group, and multiline queues. In 
the present paper we study properties of the limit distribution T and find 
some remarkable properties of it. For simplicity, most of our results are 
focused on the case when the n particles are labelled by {1,... ,n}, so we 
have a cyclic permutation. 

If we condition on the permutation n of the distribution we can in some 
special cases give an exact description of the density function g n of how the 
particles are located on the circle, see Section [3j For the reverse permutation 
wq = n... 321 the density function is the Vandermonde determinant. For a 
number of other permutations the density function is, mostly conjecturally, 
a sum of derivatives of the Vandermonde determinant. In the cases we can 
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prove, we first prove an exact formula for the discrete case with a given 
number of empty sites and then take the limit. An interesting observation is 
that the the density functions in several cases satisfy the Laplace equation. 

The probability for the particles to form a given permutation seems in 
general difficult to compute and we have no general conjecture. It can be 
computed for wo, see Theorem 13.91 In Section 0] we study the probability 
of two given labels being next to each other. There is a tantalising pattern 
for this correlation that we formulate as a general conjecture. We prove the 
conjecture in a few special cases. In Section [5] we study the probability that 
k particles adjacent to each other form a descending sequence. This is also 
proven to be the Vandermonde determinant. The proof is in the finite case 
and it proves Theorem 8.1 announced in [9]. We don’t know if these two 
different probabilities both given by the Vandermonde determinant are just 
coincidences or if there is a deeper connection. See beginning of Section 
[5] for a discussion. In Section [6] we end with some speculations about the 
similar nature of this model and the famous Razumov-Stroganov model for 
linking patterns in a disk. 

Acknowledgement We thank Andrea Sportiello for asking good ques¬ 
tions and Philippe Di Francesco for interesting discussions on the Razumov- 
Stroganov chain. 


2. Background and definitions 

2.1. Multi-type TASEP on a ring. Consider a vector m = (mi,..., m n ) 
of positive integers and a ring with N > ^ rrii sites labelled {0,..., N — 1. 
A state of the m-TASEP chain is a placement of mi particles labelled i on 
the ring 1 < i < n such that no two particles are on the same site. A 
transition in this chain is one particle being chosen uniformly at random 
and the particle tries to jump left. If the site to the left is vacant, it jumps 
there. If the site to the left has a particle with a larger label they trade 
places. If the site to the left has a particle with a smaller or equal label no 
jump will happen. So we can think of vacancies as particles labelled n + 1, 
but the notation becomes simpler by not doing so. 

Exclusion processes has been studied extensively in general and this cyclic 
m-TASEP has been considered by several authors 

Both Matrix Ansatz solutions and more combinatorial solutions have been 
suggested. Let T m (V) be the stationary distribution of the TASEP, then 
we define T as the limit when m is fixed and N tends to infinity while 
scaling the ring to have length 1. Note that we define a limit of stationary 
distributions and not the limit of the TASEP itself. We leave that as an 
interesting challenge. 

2.2. Multiline queues. We will make extensive use of multiline queues 
(MLQs), originally defined by Ferrari and Martin [T2]. We will distinguish 
between discrete MLQs and continuous MLQs. 

A discrete MLQ of type m = (mi,..., m n ) is an (n + 1) x N array, with 
mi + • • • + m, boxes in row i for 1 < i < n. Given such an array, there is a 
labelling procedure which assigns a label to each box. See Figure I!h2l for an 
example. We label the boxes row by row from top to bottom. Suppose we 


CONTINUOUS MULTI-LINE QUEUES AND THE TASEP 


3 


o o 

OOP 

o o o o 

Figure 1. A discrete multiline queue 



Figure 2. The labelling of the multiline queue in Figure [2721 



Figure 3. A continuous multiline queue 


have just labeled row i. Pick any order of the boxes in row i such that boxes 
with smaller label come before boxes with larger label. Now go through 
the boxes in this order. When considering a box labeled k, find the first 
unlabelled box in row (* + 1), going weakly to the right (cyclically) from the 
column of the box, and label that box k. When this is done, some boxes (in 
total mj+i) remain unlabelled in row i +1. Label these i +1. Thus all boxes 
in the first row are labelled 1. With a k-bully path we will mean the path of 
a label k from its starting position on row k directly down and then along 
the k + 1st row to its box with label k , and so forth all the way down to 
the bottom row. If two fc-bully-paths are arriving at the same label k it is 
not well defined which one turns downwards and which one that continues 
on the row, but it will not matter for our purposes. 

A continuous MLQ of type m is a sequence of n ’continuous’ rows with 
mi + ■ ■ ■ + nrii boxes in row i. In this case we consider the location of the 
boxes to be numbers in the continuous interval [0,1). We label the boxes 
using the same labelling procedure as for discrete MLQs (we disregard the 
set of measure zero where two or more of the positions of the boxes coincide). 
Clearly, continuous MLQs are limits of discrete MLQs. 

The motivation of these definitions is the following theorem. 
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Figure 4. The output of the process of the last row when 
the top row together with the circles in the second row are 
given. 


Theorem 2.1 (Ferrari-Martin). Let X be m-TASEP distributed word. Then 


P[X = u\ = 


n v 


n r ( n v 

lli=l VraiH-b mi) 

where n u is the number of m-MLQs whose bottom row is labelled u. 


From this theorem it follows that T is the distribution on the bottom row 
for a uniformly chosen continuous MLQ. 

2.3. The process of the last row. Theorem 12.II can easily be extended to 
the case where m n = 0, as we now explain. Consider a (mi,... , m n _i,0)- 
MLQ. By the theorem, the labelling of row n — 1 has TASEP distribution. 
Using the proof of [12], it is easy to show that row n also has TASEP 
distribution, of the same type. It follows that the probabilistic map from 
the (re — l)st row to the reth row represents another Markov chain with 
the same distribution! Here is an example of how we will use this fact. 
Consider Figure 12.31 In the top row we have sampled a word from the 
TASEP distirbution, and on the second row we have selected 4 positions for 
the boxes uniformly at random. Then we use the same labelling procedure 
as previously. The claim, then, is that the bottom row also has TASEP 
distribution. We have been deliberately vague as to whether we are talking 
about discrete or continuous MLQs here - of course, both interpretations 
are valicQ. 


3. Discrete and continuous density functions 

In this section we will assume that m = (1,.... 1), so the particles form 
a permutation. Our main interest lies in the following quantities. 

(a) For a permutation vr, the number G n (bi ,..., b n \N ) of discrete MLQs 
of length N such that the letters in the bottom row are 7Ti,..., 7r n , 
at positions b\ < ■ ■ ■ < b n . For fixed N, summing G n (bi, ..., b n \ N ) 
over all permutations n and all increasing sequences b\ < • • • < b n , 
we get Zn = (J) (• • • (J), the total number of discrete MLQs. 

(b) We get the corresponding continuous probability density function 
for 0 < q\ < ■ ■ ■ < q n < 1 as a limit, 

9n(qi,- • • ,q n ) = Jim ■q 1 \,...,[N ■ q n \\N). 

N—t oo Ztv 

The number g n can also be obtained by a finite computation by 
looking at the (!()! permutations the boxes in the MLQ can form 

^Strictly speaking, we proved it for discrete MLQs, and it is trivial to deduce from this 
that the continuous version is also valid. 
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Figure 5. Counting multiline queues with lattice paths. 

(when ordered from the left to right on a single line), or by taking 
the part of highest degree in G n appropriately normalized. 

(c) The probability p n = f 0<qi< ... <qn<1 5 tt(<?i, • • ■, q n )dqi ... dq n that the 
letters at the bottom of a random continuous MLQ form the permu¬ 
tation 7 r. 

3.1. Probabilities of the discrete chain. We are not able to give an exact 
formula for G n except for a few special permutations and conjecturally for 
some more. 

As an example one can by studying two row MLQs see that G\ 2 (b\,b 2 ',N) = 
N - b 2 + b\ and G 2 \(b \, b 2 \N) = b 2 — h. 

More generally, we have the following formula for the reverse permutation 
wq = n(n — 1)... 1. 

Proposition 3.1. For any N > n > 2 we have 
G W0 (b l ,...,b n -,N) = detiY 6 * + ^ ^) = JJ (k - b k ) 

^ ^ ' ' l<h3<n l <k<l<n d= 1 C ' 

Proof. Suppose we have a discrete MLQ whose bottom row n is labelled by 
the reverse permutation, with a particle labeled n + 1 — i at position bi, for 
1 < i < n, where b\ < ■ ■ ■ < b n . It is a direct consequence of the construction 
of MLQs that the position of the boxes on row n — 1, b[ < ■ ■ ■ < b' n _ 1 must 
be such that bi < b' x < bi + \ for 1 < i < n and they must also correspond to 
the reverse permutation (of length n— 1). It follows by induction that each 
row is labelled by a reverse permutation. 

Thus the bully paths will never touch each other and we may use the 
Lindstrom-Gessel-Viennot lemma, see e.g. [22 . Chapter 2.7]. We may say 
that each bully path starts at the beginning of the row, that is on positions 
(r, 0) for 1 < r < n (matrix notation). See Figure [5] for an illustration. 

The number of (non-cyclic) paths using only right and down steps from 
(r, 0) to ( n,bi ) is ( Setting j = n + 1 — r gives the determinant 

in the proposition. The factor n,/=i if can be taken out and using column 
operations we reduce to the standard form of the Vandermonde determinant. 

□ 
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Next we study the permutations SkW o = n ... (A; + 2)A;(A; + 1)(A; — 1)... 21, 
where the numbers k and k + 1 have switched places in wq. For 1 < k < n, 
let Ah be the matrix with entries in rows 1 < i < n — k and entries 

( b,: /- 2 ~ 2 ) rows n — k < i < n. We conjecture the following. 


Conjecture 3.2. For N > n > k > 1 and 0<bi<---<b n <N — 1 the 
number of MLQs with bottom row SfcWo is 


G Sk w o (Pi > ■ ■ ■ j bm AT) 


00 det ' 4 ‘ 


G 


WQ " 


We can prove this conjecture for k = 1,2. 


Theorem 3.3. For k = 1,2, Coniecture \3.S\ is true, that is, 


G 


SkW 0 



det Alfc 


G 


WQ ■ 


Proof. We will first consider the case k = 1. We distinguish between two 
different types of MLQs that can result in the permutation siwo, depending 
on whether the bully path for the first class particle wrapping from rightmost 
position to the leftmost. 

The first type is an MLQ q where there is no bully path wrapping. Then 
the bully paths for the first and second class particles will touch at some 
point (r, c) and after that first class particles path will be below or on the 
path of the second class particle. We could also describe this as the permu¬ 
tations on the first r — 1 rows of q is the reverse permutation and after that 
the 1 and 2 has switched places. 

To count MLQs of the first type we define an injection into sets of cer¬ 
tain non-intersecting paths V = {L n ,..., L\}. The path L\ is formed by 
concatenating the bully path of the first class particle of q from (1,0) to 
(r, c) with the bully path of second class particle from (r, c) to (n, b n ) and 
then lifting the resulting path one step upwards. So, L\ is a path from (0,0) 
to (n — l,b n ) passing (r — l,c). The path L 2 is formed by concatenating 
the second class particles bully path from (2, 0) to (r, c ) with the first class 
particles path from (r, c) to (n, b n - 1 ). For 3 < % < n, Li is the bully path for 
the ith class particle. Thus V is a set of non-intersecting lattice paths with 
starting positions <Si = {(n, 0),..., (3,0), (2,0), (0,0)} and ending positions 
A4i = {(n, bi ),..., (n, b n - 2 ), (n, b n - 1 ), (n — 1, b n )}. All such sets of paths 
are counted by the determinant of the matrix 

Sets of paths V where the vertical distance between L\ and L 2 is always 
two or more do not come from an MLQ q of the first type. To subtract 
this overcount we need to count the number of non-intersecting paths from 
S = {(n, 0),..., (2,0), (1,0)} to Ai = {(n, bi ),..., (n, 6 n _i), (n, b n )}, which 
is preciesly G Wo by Proposition 13.11 Thus the number of MLQs of the first 
type is det(-Di) - G Wo . 

The second type of MLQs q that project to the permutation si?uo has the 
bully path of the first class particle wrapping on row two before finding a box 
to label 1, see Figure[6l The bully path of the first class particle will overlap 
with the bully path of the second class particle in the beginning of row two 
and possibly more later. Apart from that, no paths are touching. We now 
desrcibe a bijection from such MLQs to non-intersecting lattice paths V = 
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{ L n ,... ,Li, Lq} with starting positions S 2 = {(n, 0),... , (2, 0), (1, 0), (0, 0)} 
and ending positions M 2 = {(n, &i),..., (n, 6 n _ 2 ), (n, b n -±), (n—1, 6 n ), (1, IV- 
1)}. We define To as a translation one step upwards of the bully path for 
the first class particle from (1,0) until it comes to (2, N — 1). The path L\ is 
a translation one step upwards of the bully path of the second class particle. 
L 2 is the part of the bully path of the first class particle starting at (2,0). 
For 3 < i < n, Li is the bully path for the ith class particle. 
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Figure 6. On top a schematic image of an MLQ projecting 
to sitco of type 2. Below the corresponding set of n + 1 
non-intersecting lattice paths as in the proof of Theorem 13.31 

By the Lindstom-Gessel-Viennot Lemma all such sets of non-intersecting 
paths are counted by the determinant of the (n + 1) x (n + l)-matrix 
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Hence the total number of MLQs is 

Gs^wq (b 1 , ..., b n ] N ) — det D 2 T det D 1 G w 0 . 
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Expanding D 2 along the bottom row gives det £2 = N ■ det A\ — det D 1 and 
the statement for A; = 1 follows. 

The case k = 2 is very similar but more complicated. This time there are 
three different types of MLQs: no bully paths are wrapping, the second class 
path is wrapping on row 3, and the third type is when the second class path 
is wrapping on row 3 and first class path is wrapping on row 2. These are 
all the cases since if the 1 is wrapping then the 2 must also be wrapping for 
the 1 to end up last in the permutation. For each type we give an injection 
to a set of tuples of non-intersecting paths, which can be counted using the 
Lindstrom-Gessel-Viennot Lemma. 

Type 1: Assume the bully paths of the second and third class par¬ 
ticle intersect for the first time at (r, c). No other bully paths touch. 
We map such MLQs injectively to V = {L n ,...,L 1 }, with starting po¬ 
sitions Si = { (77,, 0 ),... , (3,0), (1,0), (0,0)} and ending positions M\ = 
{(n, bi ),..., (n, b n - 2 ), (n — 1, b n - 1 ), (n — 1, b n )}. Let L 1 be the bully path of 
the first class particle translated one step upwards. Let L 2 be the transla¬ 
tion one step upwards of the concatenation of the second bully path from 
(2,0) to (r, c) and the third bully path from (r, c) to (n, b n -\). Let £3 be 
the concatenation fo the remaining pieces and L* is the ith bully path for 
4 < i < n. The total number of such n-tupples of non-intersecting paths is 
the determinant of the following n x n-matrix 

As in the type 1 case above we must subtract those where the vertical 
distance between L 2 and £3 is at least 2 all the time, which again is G wo . 

Type 2: Let (r, c) be the point where the bully paths for the first and 
second class particle intersect the first time. The wrapping second bully 
path and the third bully path will overlap in the beginning of row 3. We 
map such MLQs injectively to V 2 = {£«,•••, £ 1 , £o}> with starting po¬ 
sitions S 2 = {(n, 0 ),..., ( 2 , 0 ), (1, 0 ), (— 1 , 0 )} and ending positions M 2 = 
{(n, h ),..., (n, b n - 2 ), (n - 1, 1 ), (n - 1, b n ), (1, N - 1)}. 

£0 is a translation two steps upwards of the concatenation of the first 
bully path to (r, c) and the bully path of the second class particle from (r, c) 
to (3, N — 1) before it is wrapping. £1 is a translation one step upwards of 
the concatenation of the second bully path to (r, c) and the first bully path 
from (r, c) to ( n,b n ). Let L 2 be a translation one step upwards of the bully 
path of the thrid class particle and let £3 be the bully path of the second 
class particle after it has cycled, that is, from (3,0) to (n, 6 n _ 2 ). To count 
the number of MLQs of type 2 we have to subtract of the sets V 2 where the 
vertical distance between Lq and L\ is 2 or more at all times, which can 
be coutned be lowering the start and endpoints of Lq. This means that the 
number of MLQs is counted by the difference det (£ 2 ) — det(E' 2 ), where £ 2 , 
E 2 are the following (n + 1) x (n + l)-matrices 

Type 3: The MLQs where both the first and the second bully path 
wraps around can similarly be bijectively mapped to n + 2 -tupples of non¬ 
intersecting paths starting in positions £3 = {(n, 0 ),..., ( 1 , 0 ), ( 0 , 0 ), (— 1 , 0 )} 
and ending positions M 3 = {(n, b ±),..., (n, b n - 2 ), (n — 1 , b n - 1 ), (n - 1 , b n ), 
(1, N — 1), (0, N — 1)}. These are counted by the determinant of the matrix, 
Expanding the matrices E \, £ 2 , E ' 2 , £3 along the bottom rows most terms 
cancel and give the claimed result. □ 
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We will also state the following more general conjecture for commuting 
simple reflections. For a partition k = (k\ > ■ ■ ■ k r > 1), let k' denote the 
conjugate partition. For a subset S C [r], let k(£) be the partition consisting 
of the parts ki,i £ S. With k(S')' we denote the conjugate of k(S). Let As 
be the matrix with entries '~(S )^) > where ki(S)' = 0 if i > fc m inS- 

Conjecture 3.4. For N > n and k such that n > k\ > k 2 + 1 > k 3 + 2 > 
■■■ > k r + r — 1 > r — 1 and 0 < h\ < ■ ■ ■ < b n < N — 1, the number of 
MLQs is 

G Ski ... SkrW 0 (b 1 ,...,b n -N) = £ (-l)' 5 ' H (*) det As. 

SC[r] iGS ^ 

Note that A$ is the Vandermonde matrix so it specializes to Conjecture 
13.21 for r = 1. 

3.2. Probability density functions for the continuous chain. We will 

now use the above results on the discrete chain to understand the probability 
density function g n for the continuous distribution T when the particles form 
a certain permutation 7 r. 

By setting % = h/N in the formula in definition (b) and letting N —> oo 
we can from Proposition 13.11 deduce the following. 

Corollary 3.5. For any n > 2 and 0 < q\ < ■ ■ ■ < q n < 1 we have 

9w 0 (Qi, ■ ■ ■, Qn) =n\ (qi - q k ) 

1 <k<l<n 

Similarly, Conjecture 13.21 translates to the following. 

Conjecture 3.6. For any n > k > 1 and any 0 < q\ < ■ ■ ■ < q n < 1 we 
have ^ 

9 **”» = “ ') 9 “»' 

Corollary 3.7. Coniecture \3.6\ is true for k = 1 , 2. 

Example. 

#4321 = TT (Qj - qi), #4312 = ( W-1 ) #4321 

1 <i<j<A V #4 ) 

(i d 2 \ (i 9 3 

#4231 — -U-1 #4321, #3421 — TT7-U-X- 

\ 2 dq 3 dq 4 ) \6 dq 2 dq 3 dq 4 

The polynomials g w (- ■ ■) satisfy some interesting relations whose general 
form we have not been able to pin down exactly. For example, for n < 4, all 
the g w (x i,... ,x n ) satisfy Laplace’s equation 

d 2 g w 9 2 g w _ d 2 g w _ 

dxj dx\ dx\ 

It’s a classical fact [ 6 ] that any such harmonic polynomial can be expressed 
as a linear combination of partial derivatives of a Vandermonde determinant 
(the converse, that any such combination is harmonic, is immediate). Since 


1 #4321- 
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g Wo is a Vandermonde determinant, it follows that there is some linear com¬ 
bination of its partial derivatives whose value is g w for each w. In each case 
there seems to be a particularly simple one, for example: 


( d 1 d 2 \ 

- V + dil + 2 0?J 3321 

/ a i a 2 i a 3 \ 

91432 = V 94321 

/a d 1 d 2 \ 

54132 = 1 - 7j-a-^ no-b - 54321 

V oq 3 dq 4 2 dq 3 dq 4 J 

( d 1 d 2 \ 

94213 = f 1 ‘ 8^ + 2 SjJ 94321 

For n = 5, we have calculated that 15 of the 24 polynomials satisfy 
Laplace’s equation (cyclicity modded out). So there cannot be any lin¬ 
ear recursion using differentiation operators between the polynomials g w in 
general. Perhaps there is some other set of operators which coincide with 
differentiation for small n and does extend to larger n? 

The reader may note that the part of maximal degree in g u appears to be 
±g Wo , where wq = 4321, and we choose + if and only if £(wq) ~£{u) is even. 

The implication of Conjecture 13.41 to the continuous distributions is the 
following. 


Conjecture 3.8. For k such that n > k\ > k- 2 + 1 > k 3 + 2 >••• > 
k r + r — 1 > r — 1 and any 0 < q 3 < ■ ■ ■ < q n < 1, 


5sfc^ ...Sj Cr 'WQ 


1 d kr 
k r l dq n —] tr j r \ ... dq n 



Q s k-± • • - s k r _i w o * 


An example, where the conjecture is true, is n = 4, k± = 3, &2 = 1, 
/ 1 d 3 d 1 9 4 \ 

53412 — 1 - x"o— q —a -b-^ Fb—b— 02 — 54321- 

V 6 oq 20 q 3 oq 4l dq 4 6 dq 2 dq 3 d z q 4 ) 


3.3. Probability of given permutation. One obvious question to ask 
about the distribution T is the probability p n that the particles form a 
certain permutation 7 r. We can compute the exact probability of the reverse 
permutation. 


Theorem 3.9. The probability that the particles form the reverse permuta¬ 
tion is 

1 

P w 0 - i-rn-l /2fc+l\ ' 

Hfe=l 1 fc+l ) 

Proof. If the labelling of the boxes on the last row of an MLQ is w 0 then also 
the word on the row above have to be labelled by the reverse permtuation, 
as noted in the proof of Proposition 13.11 and this holds inductively for all 
rows. All the relative positions are however not fixed since rows two or more 
apart could be on either side of each other. For example for n = 3 the box 
on the top row could be to the left or right of the middle box in the last 
row. Thus giving 2 different possible relative positions of all boxes. Such 
triangular arrays of boxes have been studied previously in other contexts 
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i\j 

1 

2 

3 

4 

5 

6 

1 

0 

1/2 

1/6 

2/15 

6/55 

1/11 

2 

1/14 

0 

25/42 

2/15 

6/55 

1/11 

3 

5/42 

1/21 

0 

19/30 

6/55 

1/11 

4 

16/105 

17/210 

1/30 

0 

106/165 

1/11 

5 

68/385 

81/770 

19/330 

4/165 

0 

7/11 

6 

37/77 

41/154 

34/231 

5/66 

1/33 

0 


Table 1. Table over Cjj(n), for n = 6. 


and are called e.g. Gelfand-Tsetlin patterns, see [20]. The enumeration of 
all such patterns seems to have been done first in [23] where it is proven 
that the number of such patterns is 

na'pi+i)!' 

Think of the (”^ 1 ) boxes in the MLQ as chosen in the interval [0,1), and 
then selecting which boxes ends up on which line. Then there are a total 

( \ 2 ' ) ways to do this. Dividing the number of Gelfand Tsetlin patterns 
with the total numbers gives the formula claimed. □ 

We have computed p n for all permutations for n < 6 and we can un¬ 
fortunately not see any obivous general pattern. For instance, (i) p Wo is 
the smallest of the probabilities and p t( i is the largest, but P-k/Pwq is not 
an integer in general, and (ii) the chain is not symmetric: p«, o7r Wo 7^ Pn in 
general. 


4. Correlations 

Even though the probability p n of a given permutation ir to appear at 
stationarity seems difficult to describe in general, the correlation of two ad¬ 
jacent elements seems to exhibit very interesting patterns in the continuous 
TASEP. Let Cjj(n) = P(rr a = i,w a +1 = j) for some a. See Table [T] for the 
values of Ci t j(6 ). 

The most obvious observation is that the columns in the upper right part 
seem to be constant. To be more precise: 

Conjecture 4.1. For every i + 1 < j we have Cij(n) = n /( n ^ J ). 

From this it would also follow that c n _= (n + l)/(2 n — 1) and cp 2 = 
4/ (n + 2). 

It seems the denominator always is a product of small primes. The data 
for n < 6 suggest in fact a conjecture covering all the c t .jS. Our main 
conjecture for the correlations in the continuous TASEP on a ring is the 
following. 
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Conjecture 4.2. For n > 2, we have the following two-point correlations 
at stationarity 


Ci,j(n ) 


n 

{ n V) ’ 

n | ni 

( n t J ) (T)’ 

rt _ n 

rc(j+l) _ n(j-l) _ n 

m rr? (?) 


*/ i + 1 < j < n, 
ifi + l=j<n, 
if j <i <n, 
if j < i = n. 


So, for any j < n the most likely situation is that it is (directly) to the 
right of j — 1. If j is very close to n then this will happen more than half 
the time. For small j (1 < j < n/6 ) the second most likely thing is that j 
is to the right of n. The most unlikely of all events is that n is to the right 
of rt — 1. 


Remark 4.3. According to the conjecture, n is always a factor for the prob¬ 
abilities. It could be tempting to divide with n and say that we are interested 
in the case when w\ = i,W 2 = j ■ This is however not true. The spacing be¬ 
tween the particles is not uniform and hence it is not uniform which particle 
is first in a given [0,1) interval. 

We can prove a few of the cases of this conjecture. 


Proposition 

(!) C2,l = 

( 2 ) C12 = 

(3) c n ^ n _\ 


4.4. The following two-point correlations hold for any n > 3. 

n _ n _4_ 

( n + 1 ) (™+ 2 ) _ (n+l)(n+2) 

n | n 4 

(T) (T) “ (- +2 ) 

_ n 2 n(n- 2) n _ 3 

- fY? fr? TT) - (2n—l)(2n-3) 


Proof. For the first two statements we will use the process of the last row. 
Assume that the particles of classes 2 and 1 are positioned at t /2 and q\ 
respectively before the process of the last row. By rotation we may assume 
that q 2 < q\. The only way to obtain a 2 followed by a 1 after the process 
of the last row is to have exactly one particle in the interval [r/ 2 , qi}- That is 
easy to see since the class 1 particle will land at the first available position 
and after that the class 2 particle will do the same. Let y = 1 — (q\ — qf). We 
know by Corollary 13.51 that <721 (<?i, < 72 ) = 2(^2 ~ <?i) = 2(1 — y). If we think 
of this as the limit of the stationary distribution of TASEP it is clear that 
the particles of classes higher than 2 will not influence the relative position 
of 2 and 1. The probability that there is exactly one of the n particles in 
the interval [< 72 , <71 ] is (™)(1 — y)y n ~ 1 . We thus obtain 

I 1 4 

c 2 ,i = P(2 followed by 1) = / 2n(l - y) 2 y n ^dy = - ———— 

Jo K n + 1 -)\ n + z ) 

The computation of c\ 2 is similar. This time there are three possibilities 
to get a 1 followed by a 2; either there are no particles in the interval [q 2 , qi] 
or all the particles are in this interval or there are all particles but one in 
the interval. Summing these three integrals gives the desired formula. Note 
that the method used above could in principle (but it quickly gets more 
complicated) be extended to Cjj, for i, j < x, if we know all the density 
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functions g n for all permutations of length at most x. The reader is invited 
to try e.g. c 3j i. 

The computation of c n , n -± is more involved. We use the continuous mul¬ 
tiline queues. As discussed in Section [3] it is only the relative positions of 
the f”^ 1 ) chosen positions in the multiline queue that are important. By 
rotation we may assume that it is the leftmost box on the bottom line that 
is given class n. We number the relative position from the right and let a mj - 
be the jth box from the right on row m, so o m>m > > ■ ■ ■ > a m .i. 

The key is to note first that for the leftmost box on row n to get class n 
there must be no box ’’queuing” at the far left which means we must have 
the inequalities a n -i,j > a n ,j for each 1 < j < n — 1 and a n ^ n > a n _ 
Secondly, for the second box from the left on row n to get class n — 1 the 
leftmost box on row n — 1 must have class to — 1 and be the only box on that 
row to the left of a nn —\ . We may summarize the inequalities as follows. 


®7l—1,71 —1 

> 

tin—2,n—2 ^ 

C^n—2,n—3 ^ 

> 

O-n-2,1 

V 


V 

V 


V 

^n,n—l 

> 

^71 — 1,71—2 ^ 

A 

7 

s 

> 

1,1 



V 

v 


V 



tin,n—2 ^ 

a n , 2 > 

> 

1 


The position of the boxes on the first n — 3 rows will not matter and we 
may resort to studying the relative position of the 3n — 3 boxes on the three 
bottom rows. From the inequalities it follows directly that a U)n = 3n —3 and 
a n - i >n _i = 3n — 4. If a n>n _i = 2n + i — 1 then a n _ 2 ,j = 2n + j — 3 for all i < 
j < n —2. The remaining entries form a standard Young tableau (transpose 
the figure of the inequalities above) with columns of length n — 2, to — 2, i. Let 
SYT n — 2 ,n— 2 ,i denote the number of standard Young tableaux of this type. 
We refer to |23j for basic facts about standard Young tableaux including 
the hook-content formula from which one may deduce that SYT n _ 2 ,n~ 2 ,i = 
( 2 n l+lMlgKzi l b _ T]- ie rotation gives a factor of 3 n — 3 and we get 

/■o„ o', \pti-2 cy T onw"-2 (2«-4+t)!(«-t)(tt-i-l) 

_ \^ n 4) 2 ^t =0 ££_ -Ln—2,n-2,i _ \ n *)• 2 ^i =0 _j]_ 

c n,n—l / 3n—3 \ ~ (3n - 4)! 

\n,n—1,71—2/ v J 

n—2 ^ n 7 2 ^ 3 

as desired. The last equality was obtained using MAPLE. □ 

5. Correlation function for initial decreasing sequence 

In this section we will prove a formula for the probability that a word 
sampled from the discrete TASEP starts with a given decreasing word. By 
a remarkable coincidence it is almost the same formula as in Proposition l3.il 
Fix the length N of the ring and let m = (1,.... 1). Suppose u is picked 

N 

from the stationary distribution of the m- TASEP. We now ask, what is 
the probability that u starts with some fixed word? In general the answer 
appears to be complicated (see [9] for words of length at most 3). However 
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in the case of a word of the type x n x n -\ .. .xi where x n > x n -i > ■ ■ ■ > X 2 , 
we show that there is a simple answer to this question. 


Theorem 5.1. Suppose u is picked from the stationary distribution of the 
m-TASEP. Fix N > x n > x n -\ > ■■■ > x^ > 1. Then, the probability 
f n (x n , • • •, X 2 ) ( 7 r for ’permutation’) that for some word v, u = x n x n -\... X 2 V, 


is 



det 



71—1 


i,j =1 


The formula in Theorem 15.11 has been previously been noted by Kirone 
Mallick [7j. One way to think about this theorem is to think of the state of 
the TASEP as a permutation matrix of size N with Is in positions (i, 7 r(z)). 
Theorem 15.II below gives the probability that the first n — 1 rows have the Is 
in positions (1, x n ),..., (n — 1, £ 2 ). By cyclic invariance the same is true for 
any n — 1 consecutive rows. Proposition 13.11 on the other hand states that 
the probability is exactly the same that the first columns of the state matrix 
have Is in positions (x n , n — 1),... , (X 2 ,1))- This is because the probability 
of the position of the smallest labels does not change if we change all labels 
n ,... ,N to just n. This interesting equality does not carry over to other 
patterns in general. 


Proof. For a vector (mi,...,m n ) of positive integers with sum N, write 
Ni = (M,;_ 1 . Mf where Mq = 0 and Mi = Ylj<i m j > 1 < * < n so that 
[N] = {1,2,..., A^} is the disjoint union of the Nf s. 

Now, write f w (m n ,m n - 1 ,..., mi) (w for ’word’) for the probability that 
a TASEP distributed word of type m starts with the word n(n — 1)... 32. 

Clearly, 

(1) fw(m n ,...,m 1 )= 'Y] ••• ’Y] U(x n ,...,x 2 ). 

i n £JV„ x 2 €N 2 

Moreover, for a fixed x n ,...,X2, the value of f n {x n ,..., X2) can be com¬ 
puted from the values of all f w (m n ,..., mi) where mi + • • • + m n = N (by 
applying Mobius inversion to equation JH)). Thus, to prove the theorem, it 
is sufficient to show that equation (PQ) is satisfied (for all m) when substitut¬ 
ing the stated formula for f n . We will now make a series of manipulations 
to the right hand side of ©, after substituting. This expression is 


T-rn-l (N\ ^2 '"22 det 

llj=l \i)x n £N n x 2 €N 2 



71 — 1 


i,j =1 


Now, note that 1 only occurs in row i and use the multilinearity of the 
determinant to move each sum inside its respective row. We get 


1 


ns (?) 


det 



ns 1 (?) 


det 


(m 


(TC 
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We will use the following determinantal identity twice. It is easy to prove 
using row operations. 

Lemma 5.2. Suppose a tJ are the entries of annxn matrix such that aij = 1 
for j = 1. Then 

det {aij) 2 j =1 = det(a(j_ 1 )Q +1 ) — aiQ'+i))”^!- 
Some further manipulation of the determinant: 


det 


Mi +1 + 1 

j 


det 


Mi + 1 

j 

n —1 


n— 1 


= det 


i,j= 1 

Mi + 1 

j ~ 1 


i,j = 1 


Mi +1 + 1\ (M\ + 1 N 

3 J V 3 / / i,j=\ 

The first equality follows from row operations. The second follows from 
Lemma l5.2h note that the first column in the matrix on the right is constant). 

Let F w (m n , m n _i,..., mi) = [I'll (m) ' fw(m n , m n - 1 ,, mi) be the 
number of multi-line queues whose bottom row starts n(n — 1)... 32 and 
has type m. 

To prove equation JT]), we should prove that 


F w (m n ,m n -i,... ,mi) = 


o 


*■*■(*) 
i =1 Vi-1/ 


det 


Mi + 1 

j ~ 1 


i,j = 1 


Move the product in the numerator into the rows of the matrix and the 
product in the denominator into the columns. Simplify the resulting expres- 

( N )( M i+ 1 ) 

Mi . i /" 1 = N+ 2 -j i.M + + 2 -j) anc ^ move the prefactors out through the 

Vj—i/ 


sion 


rows and columns again. We get 
Mi + 1 


n 

i =1 
n—1 

n 

i= 1 


N + 2 — i 

Mi + 1 

N + l-i 


det 


det 


N + 2-j 
Mi + 2 — j 

N + 2-j 
Mi + 2 — j 


i,j= 1 


i,j= 1 


This matrix has a row with all ones, allowing us to convert it back to a 
(n — 1) X (n — l)-matrix again, using (the transposed version of) Lemma [5.21 
This (and (£+£>.) - (W-j.) = (^)) yields 


( 2 ) 


n—1 


n 


Mi + 1 

N + l-i 


det 


N + l-j 

Mi + 2 — j 


n—1 

i,j =1 


It remains to show that the last expression equals F w (m n ,..., mi). Con¬ 
sider a multiline queue counted by F w (m n , ..., mi). It has n—1 rows, 
indexed 1 ,... ,n — 1 (we use matrix notation; (i, j) refers to row i and col¬ 
umn j). In row 1, 2,... , (n— 1), there are mi, mi+m 2 , • • • ,N — m n particles. 
It is easy to see that the sites in the ’triangle’ (n — 1,1),... ,(n — l,n — 1); 
(n — 2 , 2 ),... , (n — 2 , n — 1); ...; ( 2 , n — 1) are filled with boxes, and that 
a box in such a site (i.j) is labeled n + 1 — j. Denote by Zj^ the distance 
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from the right end of the multiline queue of the jth box from the right in 
the (n — z)th row. That is, if the jth particle from the right in the (n — z)th 
row is at ( i, r ), we let Zij = N — r. Since the word starts with the descend¬ 
ing sequence, no bully paths may wrap around to the beginning. Thus the 
numbers Zij must form a semi-standard Young tableau (SSYT) of shape A 
where the conjugate partition is A' = M n _j — (n — * — 1) for 1 < i < n — 1. 
Moreover this is a bijection from the MLQs counted by F w (m n ,... ,m\) to 
SSYT of shape A with entries in [t], where t = N — n + 1. 

Example. 

Here, N = 13, n = 5, m = (2, 2, 2,3), = (2,4,6,9,13),* = 

9, A' = (6,4,3,2). 

The multiline queue 



5 @©©@ 5 5 ©© 5 ®@© 

corresponds to the tableau with rows 1125, 2368, 359, 57, 6, 9. 

Now we are in a position to finish our argument. 

By the definition of the Schur function, the number of SSYT of shape A 
with entries in [t] is s^(l*). Now, recall the hook-content formula and the 
Jacobi-Trudi identity, see e.g. ITT . 


Lemma 5.3. The number of SSYT of shape A and entries in [*] equals 


n = det 


rgA 


h\{r) 


K ~i + j 


= det 


t + j-1 

K - i + J 


where for a box r = (i, j) in the Ferrers diagram of A, we let c\(r) = j — i 
and h\(r) = \ + A'- — i — j + 1. 

The first two equalities are well-known, and the last is easily obtained by 
column operations. 

So, by the lemma, 


F w (m n , ... ,mi) = ]J[ 

rgA 


N - n + 1 + c\ (r) 
h\(r) 


Now let s = t Y 1 = N + l — n and = A' + 1 for 1 < z < n — 1. It is easy 
to see that our determinant in ([2]) equals (after reversing the numbering of 
rows and columns) 


71—1 

n 

Z— 1 


Mi + 1 

N + 1 -i 


det 


n—1 


s + j-1 

Ti-i + 3j J iJ=l 


which by the lemma (temporarily letting ji and s play the roles of A and 
t) equals 









CONTINUOUS MULTI-LINE QUEUES AND THE TASEP 


17 


I r Mi + 1 tt JV-n-|-2-|-c (1 (r) 

H N + 1 — i ^ h u (r) 

To prove that F w (m n ,... ,m\) equals the expression @ it thus remains 
to show that 


tt Mi + 1 -| r N — n + 2 + c^{r) ^ -p-p N — n + 1 -\- c\(r) 

H N + 1 — i hn(r ) h\(r) 

i= 1 ren ' r€A v ' 

or, in terms of Ah 


n 

rG A 


N - n + 1 + ca(c) 
h\{c) 


71 — 1 


n X'i + n — i -p-r N — n + 2 + (c) 

At -4- 1 — i 11 


\ - ! 

i =1 r£fi 


V( C ) 


A). 


This is easily checked (think of /r as the result of adding a row on top of 


6. Razumov-Stroganov 

During this work, we were struck several times by the similarity between 
our chain and the Markov chain of Razumov and Stroganov [21]. We think 
this similarity can be an inspiration to understand the (continuous) TASEP 
better, e.g. perhaps adjusting the differentiation operator in Section lT2l so 
it becomes valid for all permutations. 

A linking pattern on [2 n] is a fixed point free non-crossing involution of 
[2n], We draw linking patterns as diagrams as in Figure [7[ Let Q n be the 
set linking patterns on [2n]. For any i £ {1,2,..., 2 n} and any pattern L, 
let e^L be the pattern obtained from L by joining i with i + 1, and L(i) with 
L{i + 1). Here we take indices modulo 2 n. We can describe the image e*L 
of a generator e* acting on L by adding a tile below L - see Figured) The 
RS chain on ft n is given by at each time step applying e% where i is chosen 
uniformly at random from { 1 , 2 ,... , 2 n}. 

The generators e* satisfy the following relations. 

(A) — e^e^—re^ 

(B) e] = a 

(C) = ejei when i,j are distinct and non-adjacent modulo 2 n 

We list some similarities between the TASEP and the RS chain. 

TASEP 

• the states are permutations 

• stationarity measure (conjecturally) is largest at identity, smallest at 
reverse identity. [T4] , [I], [12J (both TASEP and continuous TASEP) 

• sum of entries of previous stationary measure equals largest compo¬ 
nent of next [13] 

• defined in terms of the NilCoxeter algebra [T2] 

• (continuous TASEP) recursions that are rooted at the reverse iden¬ 
tity. 

• (continuous TASEP) measure at reverse identity is Vandermonde 

• the k -TASEP has the same stationary distribution for all k [18] 

Raz umov- S t r oganov 
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Figure 7. A linking pattern L on [2 n] for n = 3 with L(l) = 
4, L(2) = 3, L( 5) = 6. 



Figure 8. The linking pattern L' = L, where L is from 
Figure H We have L'{ 1) = 6, L'{ 2) = 3, L'( 4) = 5. 

• the states are linking patterns 

• statationarity measure is largest at least nested linking patterns, 
smallest at most nested patterns. 

• sum of entries of previous stationary measure equals largest compo¬ 
nent of next. [21 j 

• defined in terms of the Temperley-Lieb algebra. m 

• recursions that are rooted at the most nested pattern na. 

• measure at most nested pattern is Vandermonde na- 

• the £;-RS has the same stationary distribution for all k. m 

We now explain the last point in more detail. The fc-TASEP is a gener¬ 
alisation of the TASEP where in each time step, a A;-subset of positions is 
chosen, and then ringing a TASEP bell at each position with the rule that 
for any pair of neighbouring positions, the position to the left is activated 
before the one to the right. Remarkably, the stationary measure of this 
chain is independent of k. R makes perfect sense to define a ”/c-Razumov- 
Stroganov” chain in the same way; for a Ausubset S of [2n], define eg as the 
product of all e* for i G S, taking e, before ej+i if both i and i + 1 belong 
to S. For example, for n = 4, en 4 , 7 , 8 } = e 4 e i e 8 C 7 = eieseye^ 

Let Mn be the transition matrix of the k -RS chain on The following 

theorem is an analog of the k -TASEP from [18]. It is a consequence of 
Lemma 1 in m- 

Theorem 6.1 (DiFrancesco, Zinn-Justin). The stationary distribution of 
Ain ^ is the same for all k = 1,2,..., 2n — 1. 
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